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We study the properties of the striped phases that have 
been proposed for the doped cuprate planar quantum anti- 
ferromagnets. We invoke an effective, spatially anisotropic, 
non-linear sigma model in two space dimensions. Our theo- 
retical predictions are in quantitative agreement with recent 
experiments. We focus on (i) the staggered magnetization at 
T = and (ii) the Neel temperature, as functions of doping; 
these have been measured recently in hsL2-x Stx Cu O4 with 
< a; < 0.018. Good agreement with experiment is obtained 
using parameters determined previously and independently 
for this system. These results support the proposal that the 
low doping (antiferromagnetic) phase of the cuprates has a 
striped configuration. 
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Our understanding of the undoped planar cuprates re- 
lated to high temperature superconductors has been in- 
formed largely by the insights of Chakravarty, Halperin 
and Nelson |l) (hereafter, CHN). It is now widely ac- 
cepted that the antiferromagnetic phase of these materi- 
als can be well described by an effective non-linear sigma 
model. However, the situation is not so clear when we 
turn to the behavior of doped systems, as the antiferro- 
magnetic moves toward a superconducting instability. It 
has long been recognized ||,^ that there is a tendency for 
the holes to phase separate in these strongly correlated 
insulators. Experimental observations Q support this 
scenario for some of the La-cuprates at high concentra- 
tions of dopants. In addition, recent nuclear quadrupole 
resonance (NQR) and muon spin resonance (/iSR) exper- 
iments Q in La2_2;Sr2;Cu04, with < a; < 0.018, have 
been interpreted within a picture where holes are seg- 
regated into a set of parallel stripes. It is the primary 
goal of this Letter, in fact, to do a more careful theoret- 
ical analysis of the consequences of that model and to 
compare the predictions to those experiments. 

With the xy coordinate axes chosen as the crystal axes 
in the Cu02 plane, we take an array of equivalent uni- 
formly spaced stripes parallel to the y-axis. In this sim- 
plest of models we neglect both static and dynamic fluc- 
tuations (and domain formation). The situation is not 
very different if we consider diagonal stripes such as the 
ones found in La2Ni04.i25 and Lai.8Sro.2Ni04 But 
here we will treat only the horizontal stripe geometry 
suggested for 1^3,2- x^TxCuOi by the experimental data 
. In the regions between the stripes the Cu spin corre- 
lations are presumed to remain antiferromagnetic. 



As in the CHN analysis |^] of the pure case, we will 
argue that the underlying symmetries suggest an appro- 
priate nonlinear sigma model to describe the long wave- 
length behavior which determines the phase diagram. We 
start from a microscopic model with the spins interacting 
via Heisenberg antiferromagnetic exchange in the regions 
between stripes, and by a weaker such exchange across 
each stripe. Conceptually, we use real space renormaliza- 
tion in the x-direction to integrate out the short distance 
dynamics, starting with blocks the size of the region be- 
tween stripes. This is continued to a scale of the correla- 
tion length, much larger throughout the region of interest 
to us than the inter-stripe spacing. The expected lead- 
ing result for the low frequency dynamics in which we 
are interested is a spatially anisotropic nearest neighbor 
Heisenberg exchange model, 

^ ^ "^y S{n,m) ■ S{n,m') 

n,<m,m' > 

+ Jx S{n,m) ■ S{n' ,m) (1) 

m,<n,n'> 



where n and m label the sites in the x and y directions, 
respectively. Jx and Jy are the effective spin exchange 
strengths in the two directions. In general they depend 
on the doping x, but the exact form of that dependence 
relies, of course, on microscopic details. 

The destruction of long range order at zero tempera- 
ture and the suppression of the staggered magnetization 
with increased hole doping then come about from the di- 
mensional crossover from two to one dimensional behav- 
ior with the increasing anisotropy of the system (in our 
effective cr-model the critical coupling constant gc, which 
separates ordered from disordered behavior at T = 0, 
decreases with doping). At the same time, the renormal- 
ized spin stiffness decreases as the density of holes grows. 
As gdx) decreases it equals at some critical concentra- 
tion Xc the coupling constant of the undoped system and 
antiferromagnetism disappears. 

To simplify the study of the long wavelength physics it 
is convenient to take the continuum limit of (H). We use 
the usual spin coherent state representation |p| and we 
find that the effective action (in euclidean time) in the 
partition function can be written as 
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where n is a unit vector and a is the sub-lattice constant. 
This gives spin wave velocities Cy = 2S'^a'^Jy{Jx + Jy)/'h^ 
and = 2S'^a? Jx{Jx + Jy) /Tt-^ ^ which agree with a simple 
(noninteracting) spin wave calculation based on (|l|). This 
action was studied numerically in a different context a 
few years ago [||. 

It is useful to rewrite (^) more symmetrically by 
a dimensionless rescaling of the variables: x' = 
{Jy/ Jx^^^xA, y' = [Jx/ JyY^'^yA (A is a momentum cut- 
off), and t' ~ 'j2{Jx + Jy)yJ JxJySar /fi. Then the ef- 



fective action (g) becomes 



(250) 
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where ^ takes the values x',?/',t', qq — HcqA/p'^ ~ 
\2{Jx + Jy^l \fJxJy\ {aA)/S is the bare couphng con- 
stant, Co = [2{Jx + Jy)yJ JxJyY^'^{aS) /h the spin wave 
velocity and p° = ^ JxJyS^ the classical spin stiffness of 
the rescaled model. The original anisotropy is now hid- 
den in the limits. We started with a problem with a finite 
bandwidth, a lower bound on length which requires us to 
impose a cutoff in the original continuum formulation. 
The change of variables introduces an anisotropy in the 
cutoffs. 

The (T-model action, and the spin correlations it im- 
plies, can be studied in the large N limit {N is the number 
of components of fi) , where a saddle point approximation 
becomes exact [^,^. The staggered spin-spin correla- 
tion function is given by S(fc, cjn) = go/{k'^ + + 
where m is defined by the self-consistent condition, 
Sfc '^{^1 ^n] ITT-) = 1, and sets the scale for the corre- 
lations in the system. The self-consistency equation can 
also be written, after the sum over Matsubara frequen- 
cies, as 
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d'^k coth(Vfc2 + m2/3coA/2) 



{2n) 



VP 
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(4) 



Formally, this integral has a logarithmic ultraviolet di- 
vergence. If we choose an isotropic cutoff n/Aa for 
the Fourier transform of the original problem then 
\kx\ < 7ra^/'*/(Aa) and \ky\ < 7ra~^/'^/(Aa), where 



Jx/ Jy7 



(5) 



is the anisotropy parameter (as usual, we take Aa — 2y^ 
to preserve the area of the magnetic Brillouin zone). 
By the definition above, the coupling constant 5o(ck) = 
5o(l)\/(l + oi)/{2y/a). In what follows we will choose 
(70(1) renormalized to give the interacting spin wave ve- 
locity of Oguchi | p2[ , rather than the simple spin wave 
value above. Because the effective exchange Jx in the x- 
direction is weakened by the stripes, we have < a < 1. 

At high anisotropy {a 0) the spin chains become 
disconnected and we find one dimensional behavior. In 



this limit m » a^^^, and we find at zero temperature 
from d), m{a) « ia-^/^/^/^ex^p\-ATr^/^/{a^/^goia))] , 
which is just the Haldane gap [ [l3| . We came to this re- 
sult because the non-linear cr-model without a topological 
term describes the behavior of an integer spin chain . 
In the opposite limit, m « a^l^ , the system is fully two 
dimensional and Eq. (^) at zero temperature leads to 
K, 4:V2TTa^/^/VT+a[l/gc{a) - l/go{l)] where 



g^{a) = Stt^/ Va/(1 + a) {In {Va + Vl + a) 



+ \/^ln[(l + + a)/\/a]}" 



(6) 



is the critical coupling constant of the theory. It is easy to 
see that gc decreases monotonically to zero as the system 
becomes more anisotropic and the increasing quantum 
fluctuations again make magnetic order less stable. Di- 
mensional crossover occurs when the correlation length 
in the x direction becomes of the order of the lattice 
spacing, that is, m{at) ~ a^^, and a numerical estimate 
gives at ~ 0.001. 

To understand the corrections to the classical limit we 
employ a renormalization group (RG) calculation up to 
second order in the coupling constant. Proceeding as in 
[|l|, we find within a one-loop approximation the renor- 
malized spin stiffness 



p,(a)«p°(a)(l-5o(l)/5c(a)). 



(7) 



reduced from its classical value for fixed a. As a 
first consequence we see that the stiffness will vanish 
for some critical value, a — ac, which can be calcu- 
lated numerically. Oguchi's spin wave theory [l|jl^ gives 
.go(l) ~ 9.536. Then (0) gives ac ~ 0.047, much greater 
than at, so the transition from classical to quantum 
regimes (ordered to disordered ground state) occurs well 
before the potential crossover from two to one dimen- 
sional behavior [|l^ . 

We turn to the finite temperature behavior. Again, 
we can use the inequality m a^l^ in estimating the 
correlation length ^ cx 1/m from equation (^). However, 
we know from studies of the undoped system that the one 
loop calculation does not give a very good result there for 
the prefactor of the temperature dependent exponential 
expression for ^, and the same will surely be true for 
a < 1 . Here we use an interpolation formula between the 
exact result of Hasenfratz and Niedermayer ||l^ for the 
nonlinear sigma model, which is valid close to the ordered 
phase, and the result for the renormalized critical regime 
where i oc T"! |l|: 



i{T,a) 



hco 



„27rp,(a)/fcBT 
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In general the staggered magnetization depends on the 
short, as well as on the long wavelength physics of the 
problem. As in [Q] , instead of using the non- linear sigma 
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model directly we observe that the equal time staggered 
spin-spin correlation function at large distances becomes 
proportional to the square of the staggered magnetiza- 
tion, lim2;_^oo^j^_+oo 2(a;,?/) = {Ms/M^f , where Mq is the 
staggered magnetization for the classical problem (the 
fully aligned Neel state, without quantum fluctuations). 
At zero temperature, in the ordered phase, the correla- 
tion length is infinite. The only scale left is the Josephson 
length, = Hcq/ps- On one hand the asymptotic (large 
distance) limit of the correlation function has been es- 
tablished at this scale. On the other hand this is the 
maximum scale at which the long wavelength dominated 
theory we have developed can be trusted. Therefore, we 
equate the correlation function at this (scaled) length to 
the square of the relative zero point magnetization and 
using the RG result (0), we find 

Ms{a) _ l l-goil)/gc{a) , . 

Ms{l) ^ l-go{l)/gc{l)' 

As expected, Ms{ac) — 0, and close to ac we find 

1/2 

Ms{a) w {a/ac — 1) which has the mean field expo- 
nent of 1/2. 

To have true long range order at T > we need to 
invoke the weak coupling between planes, J±. Because 
J±/J « 5 X 10~^ H^] is so small, long range correla- 
tions have built up in the plane well above the ordering 
temperature T/y. Spin fluctuations then involve large cor- 
related regions, and it is a good approximation to treat 
Jj_ within a mean field theory. To determine the depen- 
dence of the Neel temperature on anisotropy one needs to 
know the connection between the classical spin stiffness 
and a. Here we assume as a simple model that the clas- 
sical stiffness is approximately independent of a over the 
region of interest, equal to the undoped value p^{a — 1), 

so ^JJxJy — J • 

Magnetic order is destroyed by thermal fluctuations, 
at a temperature T/y, when the energy ksT^ becomes 
sufficient to flip the spins in a region of linear dimension 
of the order of the coherence length ^. Since the number 
of spins in this region is proportional to (^/a)^, and the 
relative staggered magnetization in the region is given by 
Ms/Mq, we estimate 

kBTN[a)Kjj_l —I . (10) 

This expression has been used previously |l[] to estimate 
J±/kB ~ 0.01 K from the experimental Tjv of the pure 
material. Since Ms/Mq < 1, this gives (^/a) > 10 for 
Tjv > 1 K, suggesting that this mean field theory is rea- 
sonable for T/v greater than a few kelvin. We now sub- 
stitute the expression (|^) for ^ in the quantum critical 
regime into (|l^) to find TA,(a) « {a/ac - l)^^^ defining 
a critical exponent of 1/3. 



In order to compare our results with the experiments 
we need the relationship between the anisotropy param- 
eter (H) and the doping concentration x. Explicit calcu- 
lation of that is complicated, depending intrinsically on 
the chemistry of the material; it goes beyond our previous 
analysis. Instead, we propose a simple parametrization, 

a(a;) = e-^/^«, (11) 

with xq, the single (non- universal) parameter character- 
izing the behavior, dependent in an unkown way on the 
microscopic details of the system. We fix the value of 
Xq from experiment. We stress that this is the only free 
parameter in our theory. All others are obtained inde- 
pendently for the undoped material, and they are well 
known. 

Experimentally the staggered magnetization and the 
Neel temperature vanish at a critical doping Xc ~ 0.02 — 
0.023 j5|. This value of Xc must correspond in our model 
to the point at which the spin wave stiffness vanishes, 
that is, to the critical value ac ~ 0.047. Then ( pT| ) gives 
Xq w 0.007. Depending on the density of holes in the 
stripes, this corresponds to a stripe spacing of order 10 
lattice constants or so, a reasonable scale on which to ex- 
pect substantial changes to occur. Substituting ( pl| ) into 
(^) , we find the staggered magnetization as a function of 
doping. In Fig. (^ we show the theoretical prediction 
(full line) and the experimental values (dots) from ref. 

In view of the approximations used in the theory 
above, the agreement is impressive. 

M (x)/M (0) 
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FIG. 1. Staggered magnetization (normalized relative to 
the undoped case) as a function of doping (line) and experi- 
ment (dots). 



3 



The calculation of the Neel temperature from ( |lO| ) 
brings in new parameters, including the interplanar ex- 
change J±_ , but they are all set by those of the undoped 
material. We need no additional assumptions; we simply 
use the same parameters taken by CHN to fit exper- 
iment: Psil)/kB ~ 187 K and hcQ{l) « 7.7 x 10^ A° K. 
Using equations (||) and (|l^) we obtain the full curve in 
Fig. ^. The agreement between theory and experiment 
is very good except for the smallest value of Tjv, and 
can be improved by small adjustments of the parameters 
p8[ . Our conclusion is, therefore, that our theory is fully 
consistent with the recent experiments on ha,2-x Sr^, Cu 
O4. 



T (x)/T (0) 




X 



0.005 0.01 0.015 0.02 



FIG. 2. Neel temperature (normalized relative to the un- 
doped case) as a function of doping. Theory (continuous line) 
and experiment (dots). 

In summary, we have developed an effective field theory 
for magnetic correlations in doped antiferromagnets that 
are anisotropic due to the stripe correlations expected 
from microscopic phase segregation of holes. We suggest 
that the long wavelength spin fluctuations are described 
by a spatially anisotropic non-linear sigma model. The 
experimental parameter, doping concentration x, is re- 
flected in the theory by the magnetic exchange anisotropy 
a, taken to depend exponentially on x (introducing the 
only free parameter of the theory, the decay rate xo)- We 
further assume that the classical spin stiffness is not af- 
fected by doping. We believe that the good agreement 
between this theory and experiment gives new support 
to the picture of stripe correlations of holes in La2-x Sr^^ 
Cu O4 for < a; < 0.018. 
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